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Abstract. We discuss a general method by which a higher order difference equation on a group 
is transformed into an equivalent triangular system of two difference equations of lower orders. 
This breakdown into lower order equations is based on the existence of a semiconjugate relation 
between the unfolding map of the difference equation and a lower dimensional mapping that 
unfolds a lower order difference equation. Substantial classes of difference equations are shown to 
possess this property and for these types of equations reductions of order are obtained. In some 
cases a complete semiconjugate factorization into a triangular system of first order equations is 
possible. 

1 Introduction 

There are a number of known methods by which a difference equation can be transformed into 
equations with lower orders. These methods are almost universally adapted from the theory of 
differential equations and can be divided into two major categories. One category uses symmetries 
in the solutions of a difference equation to obtain coordinate transformations that lead to a reduction 
in order; see, e.g., [3], [10], [13], [H]. The second category consists of methods (e.g., operator 
methods) that rely on the algebraic and analytical properties of the difference equation itself; see, 
e.g., [7]> [II]- For historical reasons the first of these categories has been developed considerably 
more than the second. 

In this article we introduce a general new method that falls into the second category. This 
method uses the existence of a semiconjugate relation to break down a given difference equation 
into an equivalent pair of lower order equations whose orders add up to the order of the original 
equation. A salient feature of the equivalent system of two equations is its "triangular" nature; see 
PQ for a formal definition as well as a study of the periodic structure of solutions. Also see |26j for 
additional information and a list of references. Such a system is uncoupled in the sense that one of 
the two equations is independent of the other. 

In some cases (e.g., linear equations) we use the semiconjugate relation repeatedly to transform 
a higher order difference equation into an equivalent system of first order equations with triangular 
uncoupling. In the specific case of linear equations, the system of first order equations explicitly 
reveals all of the eigenvalues of the higher order equation. 
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Semiconjugate-based reduction in order uses only patterns or symmetries inherent in the form 
of the difference equation itself. For this reason, we do not limit ourselves to equations defined 
on real numbers. Allowing more general types of objects not only improves the applicability of 
the method to areas where discrete modeling naturally occurs (e.g., dynamics on networks, see for 
instance, |15j ) but it is actually helpful in identifying form symmetries even within the context of 
real numbers where more than one algebraic operation is defined. After defining the basic idea of 
semiconjugate factorization and laying some groundwork we discuss applications of the method to 
reduction of order in several classes of difference equations. 

2 Difference equations and semiconjugacy 

Difference equations of order greater than one that are of the following type 

Xn+l = fn(%m Xfi—lj • • • , 3^n— fc) (1) 

determine the forward evolution of a variable x n in discrete time. As is customary, the time index 
or the independent variable n is integer- valued, i.e., n € Z. The number A; is a fixed positive integer 
and k + 1 represents the order of the difference equation ([I]). If the underlying space of variables 
x n is labled G then f n : G k+1 — > G is a given function for each n > 1. If f n = / does not explicitly 
depend on n then dU) is said to be autonomous. A solution of Eq.([T]) is any function x : Z — > G 
that satisfies the difference equation. A forward solution is a sequence {x n }^L_ k that is recursively 
generated by ([1]) from a set of k + 1 initial values xq,x-i, . . . , £ G. Forward solutions have 
traditionally been of greater interest in discrete models that are based on Eq.([T]). 

If the functions f n do not have any singularities in G for all n then clearly the existence of 
forward solutions is guaranteed from any given set of initial values. Otherwise, sets of initial values 
that upon a finite number of iterations lead to a singularity of f n for some n must either be specified 
or avoided. In this paper we do not discuss the issue of singularities broadly but deal with them 
as they occur in various examples. 

We assume that the underlying space G of variables x n is endowed with an algebraic group 
structure; for a study of semiconjugacy in a topological context see [23] or [25]. Often there are 
multiple group operations defined on a given set G. This feature turns out to be quite useful. 
Compatible topologies (making group operation(s) continuous) may exist on G and often occur 
in applications. Topological or metric concepts may be used explicitly in defining certain form 
symmetries and are required in discussing the asymptotic behaviors of solutions with unknown 
closed forms in infinite sets. 

We "unfold" each f n in the usual way using the functions F n : G k+1 — > G k+1 defined as 

F n (u , ...,u k ) = [f n (uo, . . .,u k ),uo, . . . ,u k -i], Uj e G for j = 0, 1, . . . ,k. (2) 
The unfoldings (or associated vector maps) F n determine the equation 

(2/Q,n+l; Vl,n+1, ■ ■ ■ i Vk,n+l) = Fn{y0,n, Ul,n, • • • , 2/fc,n) 
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in G k+1 or equivalently, the system of equations 

Vj,n+l = 



fn(yo,n,yi,n, ■ ■ ■ ,Vk,n), if J = 

yj-i, n , for j = l,...,k 



In this context we refer to (yo,n> yi,m ■ ■ ■ 1 2/fc,n) as the state of the system at time n and to G k+1 
as the state space of the system, or of 

Definition 1 Lei S and M be arbitrary nonempty sets and let F, <3? be self maps of S and M , 
respectively. If there is a surjective ( onto ) mapping H : S — > M such that 

H o F = $ o H (3) 

then we say that the mapping F is semiconjugate to & and write F > We refer to $ as a 
semiconjugate (SC) factor of F. The function H may be called a link map. If H is a bijection 
(one to one and onto) then we call F and <!> conjugates and alternatively write 

$ = H o F o H~ l . (4) 

We use the notation F ~ $ when F and $ are conjugates. Straightforward arguments show 
that ~ is an equivalence relation for self maps and that > is a transitive relation. 

As our first example shows, the fact that the link map H is not injective or one-to-one is 
essential for reduction of order. Such a non-conjugate map often exists where conjugacy does not. 

Example 2 Consider the following self maps of the plane M?: 

F\(u, v) = [u 2 + v 2 , 2uv], F2(u, v) = [u 2 — v 2 , 2uv]. 

Note that u 2 + v 2 + 2uv = (u + v) 2 so let H\(u, v) = u + v to obtain 

H\(Fi(u, v)) = {u + vf = [H^v)} 2 . 

It follows that Fi is semiconjugate to = t 2 on M. For Fi define H2(u,v) = (u 2 + v 2 ) 1 ^ 2 to get 



H 2 (F 2 {u, v)) = ^ {u 2 - v 2 ) 2 + 4u 2 v 2 = u 2 + v 2 = [H 2 (u, v)} 2 . 

Therefore, F 2 is also semiconjugate to $(f) = t 2 though this time on [0, oo) since H 2 (u,v) > for 
all (u,v) G M 2 . 

In the next example we determine all possible semiconjugates for a difference equation on a 
finite field. 
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Example 3 Consider the following autonomous difference equation of order 2 on Z2 = {0, 1} 

x n+ x = (x n + l)x n -i, xq, x-i € {0, 1} 

Note that Z2 is a field under addition modulo 2 and ordinary multiplication of integers. The 
unfolding of the defining map f{u,v) = (u + l)v is given by F(u,v) = [(u + l)v,u}. In order to 
determine the semiconjugates of F on Z| first we list the four possible self maps 0/Z2, i.e., 

0o = 0, (pi = 1 (constant functions) 
t(t) =t (the identity function) 

4>(t) = t + i. 

Next, there are 16 link maps H : Zg — > Z2 of which the 14 non-constant ones are surjective. We 
list these 14 maps succinctly as follows: 

#1 = (0,0, 0,1) # 2 = (0,0,1,0) # 3 = (0,0,1,1) #4 = (0, 1,0,0) 

# 5 = (0,1, 0,1) ff 6 = (0,1,1,0) #7 = (0,1, 1,1) #8 = (1,0,0,0) 

#9 = (1,0, 0,1) #10 = (1,0,1,0) i?u = (1,0, 1,1) #12 = (1,1,0,0) 

#13 = (1,1, 0,1) #14 = (1,1,1,0) 

In the above list, for instance #14 is defined by the rule 

#i 4 : (0, 0) -> 1, (0, 1) - 1, (1, 0) - 1, (1, 1) - 

and can be written in algebraic form as #14 (u, i>) = m> + 1. T/ie oi/ier Zmfc maps are defined in the 
same way. Now writing F as 

F: (0,0) (0,0), (0,1) -(1,0), (1,0) — (0,1), (1,1) ^ (0,1) 

and calculating the compositions Hj o F for j = 1, ... ,14 the following results are obtained: 

H x oF = (0,0,0,0) = O °#l, #14 o F = (1, 1, 1, 1) = 0i o #14 
H 7 oF = (0, 1, 1, 1) = #7 o t , # 8 o F = (1, 0, 0, 0) = # 8 o t . 

Therefore, F is semiconjugate to each o/0o,0i and t 6at not to (Trai also see Example^ below). 

It is worth noting that the maps #7 and #§ are invariants (see, e.g., [2], [9], [E]) because 
the factor map is the identity function (l here). It is easy to check that in algebraic form the two 
invariants can be written as 

Hj(u, v) = uv + u + v, #s(n, v) = uv + u + v + 1 = (u +!)(« + 1). 
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3 Semiconjugate factorization 



In this section we define the concepts needed to present our results in subsequent sections. 
3.1 Form symmetries, factors and cofactors 

To make the transition from semiconjugate maps to reduction of order in difference equations, let 
A; be a non-negative integer and G a nonempty set. Let {F n } be a family of functions F n : G k+1 — > 
G k+1 where F n = [/i )Tl , . . . , /fc+i, n ] and fj, n '■ G — > G are the component functions of F n for all j 
and all n. 

Let m be an integer, 1 < m < k + 1 and assume that each F n is semiconjugate to a map 
<]?„ : G m — * G m . Let H : G k+1 —> G m be the link map such that for every n, 

HoF n = $ n oH. (5) 

Suppose that 



H(uq, ...,u k ) = [h\{uo, . . .,u k ), h m (u , . . .,u k )} 

(ti , ■ ■ ■ , tm), • • • , (^m,n 

(ti, . . . , tfn)] 

where hj : G k+1 — > G and ^j jn : G — > G are the corresponding component functions for j = 
1,2, ... ,m. Then identity ([5]) is equivalent to the system of equations 

hj(fl, n (uO, ■ ■ ■ ,Uk), ■ ■ ■ ,fk+l,n(uQ, ■ ■ ■ ,«fc)) = 

0j,n{hi{uo,...,Uk),...,h m (uo,.--,Uk)), j = 1,2, . . . ,m. (6) 

If the functions fj jTl are given then ([6]) is a system of functional equations whose solutions hj , <f>j >n 
give the maps H and & n . Note that if m < k + 1 then the functions on G m define a system 
with lower dimension than that defined by the functions F n on G k+l . We now proceed to find a 
solution set for ([6]). 

Let G be a group and denote its operation by *. If F n is the unfolding of the function f n as in 
([2]) and F n satisfies © then it need not follow that the maps are of scalar type similar to F n . 
To ensure that each & n is also of scalar type we define the first component function of H as 

hi(u , ...,Uk) = u * h(ui, ...,u k ) (7) 

where h : G k — > G is a function to be determined. 
Using (jTj) in the first equation in ([6]) gives 

f n (u , ...,Uk)* h(u , u k -i) = (8) 
g n (uo * h(m, . . . ,u k ), h 2 (u , . . . ,u k ), ■ ■ ■ ,h m (u , . . . ,u k )) 
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where the abbreviation g n is used for 0i jTl : G m — > G. 

Eq.([8|) is a functional equation in which the functions h, hj,g n may be determined in terms of 
the given functions f n . Our aim is to extract a scalar equation of order m such as 

tn+l — 9n(tni ■ • • ; tn— m+l) (9) 

from ([8j) in such a way that the maps $ n will be of scalar type. The basic framework for carrying 
out this process is already in place; let {x n } be a solution of Eq.([T]) and define 

tn — %n * h(X n — l , ■ ■ ■ , X n — fc) . 

Then the left hand side of © is 

* h(x n , . . . , x n _ fc+ i) = t n+ x, (10) 

which gives the initial part of the difference equation Q. In order that the right hand side of (|8|) 
coincide with that in ([9]) it is necessary to define the component functions in © as 

hj(x n , . . . 1 X n _ k ) — tn—j+l — ^n— j+l * h\X n —j) . . . , X n _fc_j_(_x), (H) 

for j = 2, . . . , to. 
Since the left hand side of (fTTj) does not depend on the terms 

•"n—k—li • • • i •En—k—j+1 

it follows that the function hj must be constant in its last few coordinates. Since h does not depend 
on j, the number of its constant coordinates is found from the last function h m . Specifically, we 
have 

hm{x rl . • • • , — * h(x n — m , . . . , X n — k , . . . , X n — k ~ m +\) (1^) 

fc— m+l variables m— 1 terms with h constant 

The preceding condition leads to the necessary restrictions on h and every hj for a consistent 
derivation of © from (jSJ). Therefore, (|12p is a consistency condition. Now from (jlip and (|12p we 
obtain for (uo> • • • j u k) G 

frj(«o, ■ ■ ■ ,Wfc) = * h(uj' u j+i ■ ■ ■ >v>j+k-m), j = l,...,m. (13) 
With this definition for the coordinate functions hj we rewrite ([5]) for later reference as 

f n (u , ...,u k )* h(uo, ... , u k _ m ) = 
g n (u * h(m, . . . ,Uk- m +i),ui * h(u 2 , . . . ,u k - m+ 2), ■ ■ ■ ,u m -i * h( u m, ■ ■ ■ ,«fc)) (14) 

We now define a natural form symmetry concept for the recursive difference equation ([1]). 
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Definition 4 The function H = [hi, . . . , h m ] whose components hj are defined by i!3\) is a (re- 
cursive) form symmetry for Eq.^. Since the range of H has a lower dimension m than the 
dimension k + 1 of its domain, we say that H is an order-reducing form symmetry. 

For each solution {x n } of (pQ) according to ©, (fTUj) and (fT3|) the sequence 

tn — x n * h[X n —l, . . . , £ n _fc_|_ m _i ) 

is a solution of ([9]) with initial values 

t-j = x-j * h(x-j-i, . . . ,x-j- k+m -i), j = 0,...,m-l. (15) 
Thus, the following pair of lower order equations are satisfied: 

), (16) 

x n+l = tn+l * [h( x n, ■ ■ ■ > x n-k+m)] 1 (17) 



where —1 denotes group inversion in G. The pair of equations (j!6|) and (|17j) is uncoupled in the 
sense that (|16jl is independent of pT|) . Such a pair forms a triangular system as defined in [1] and 
[261. 



Definition 5 Eq. hlb}) is a factor of Eq. ^ since it is derived from the semiconjuyate factor ^ri. 
Eq. that links the factor to the original equation is a cof actor of Eq.^j. We refer to the 
system of equations U6\) and fiTp as a semiconjugate (SC) factorization of Eq.^j. Note that 
orders m and k + 1 — m of [T6\) and (TF?|) respectively, add up to the order of (QP- 



3.2 The factorization theorem 

We are now ready for the following fundamental result, i.e., the equivalence of Eq.([T]) with the 
system of equations (fTB]) and (JTj 



Theorem 6 Let k > 1, 1 < m < k and suppose that there are functions h : Q k ~ m + 1 —> Q an d 
g n : G m — > G that satisfy equations and 
(a) With the order-reducing form symmetry 

H{u , . . . , Uk) = [uq * h(ui, ... , Uk+l-m), ■■■ , Um-1 * h{u m , ... , U k )} (18) 



Eq. (Op is equivalent to the SC factorization consisting of the triangular system of equations (E 
and [T7\ ). 

(b) The function H : G k+l —> G m defined by [TS\) is surjective. 

(c) For each n, the SC factor map <& n : G m — > G m in (0) is the unfolding of Eq. U6\) . In 
particular, each $ n is of scalar type. 
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Proof, (a) To show that the SC factorization system consisting of equations (|16p and (|17p is equiv- 
alent to Eq.([T|) we show that: (i) each solution {x n } of ([1]) uniquely generates a solution of (fT6|) and 
([TTj) and conversely (ii) each solution {(t n ,y n )} of the system (fTU|) and (fTTj) correseponds uniquely 
to a solution {x n } of (H|). To establish (i), let {x n } be the unique solution of (P) corresponding to 
a given set of initial values xq, . . . x_k £ G. Define the sequence 

tn = x n * h[X n —\, . . . , X n _] i j rm _i) (19) 

for n > —m + 1. Then for each n > using (|14p 



^n+i fn\p^ni • • • j •^n—k) 

— 9n(p^n * h(X n —\, • • • , ^n— fe+m— l)i • • • > m+1 * h(x n — m , . . . , In-t))* 

[/l(x n , . . . , X n _^4. m )] 

— 9n(tni ■ ■ ■ j ^n— m+l) * [^(■^n: • • • ; ^n— fc+m)] 

Therefore, g n (t n , . . . ,t n - m +i) = x n+ \ * h(x n , . . . ,x n -k+ m ) = t n +i so that is the unique 
solution of the factor equation (|16j) with initial values 

t_j = x_j * h(x-j-i, . . . , x_j_fc +m _i), j = 0, . . . , m — 1. 

Further, by (fl"9j) for n > we have x n +i = t n +i * [h(x n , . . . , x n _fc+m)] 1 so that {x n } is the 
unique solution of the cofactor equation (fTT|) with initial values y_j = x_j for i = 0, 1, . . . , k — m 
and i n as obtained above. 

To establish (ii), let {(t n ,y n )} be a solution of the factor-cofactor system with initial values 

*0j • • • ; t-m+l,y-m, ■ ■ ■ V~k ^ G. 
Note that these numbers determine y_ m +i, ■■■ ,yo through the cofactor equation 

y-j = *-i * [Wv-j-ii ■ ■ ■ iV-j-i-k+mfr 1 , i = o,...,m-i. (20) 

Now using (|14p for n > we obtain 

Vn+l = tn+l * [h(y n , ■ ■ ■ , Vn-k+m)]' 1 

= 9n{tn, ■ • • ; tn-m+l) * [h(y n , ■ ■ ■ j 2/n— fe+m)] 

= #n(ZM * h(yn-l, ■ ■ .,y n -k+m~l), ■ ■ ■ , Vn-m+1 * h(y n -m, ■ ■ ■ ,2/n-fc))* 

[h(y„, . . . ,y n _ fc+m )] _1 

= fniyni ■ ■ ■ i Vn—k) 

Thus {y n } is the unique solution of Eq.(H|) that is generated by the initial values (|20|) and 
y- m , ■ ■ ■ y~k- This completes the proof of (a). 
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(b) Choose an arbitrary point [v± , . . . , v m ] G G m and set 

Um-l =v m * h(u m ,u m+ i UkY 1 
where u m = u m +i = . . . Uk = u where u is a fixed element of G, e.g., the identity. Then 

v m = n m _i * h(u,u.. . ,u) 

= Um-l * h(u m ,U m+ i ...,Uk) 

= h m (u , ... ,u k ) 

= h m (u , . . .,u m -2,v m * h(u,u. . . ,n)~ 1 ,n. . . ,u). 
for any choice of points uo, . . . , u m -2 € G. Similarly, define 

U m -2 = fm-1 * h(u m -l,U m . .., M fe _l) _1 

so as to get 

f m _i = n m _ 2 * h(u m -x, u m . . . ,Ufc_i) 
= hm-i(no, ... ,u k ) 

= hm-i(uQ, . . . ,U m _ 3 ,V m -i * h(u m -x,u.. . ,u)~ 1 ,ti m _i,W. . . ,u) 

for any choice of uq, . . . , u m _3 G G. Continuing in this way, induction leads to selection of u m -i, ■ ■ ■ ,uo 
such that 

vj = hj(u , . . .,u m -i,u. ..,u), j = 1, . . . ,m. 

Therefore, H(uq, . . . , u m _i, u . . . , u) = [vi, . . . , v m ] and it follows that H is onto G m . 

(c) It is necessary to prove that each coordinate function <ftj tn is the projection into coordinate 
j — 1 for j > 1. Suppose that the maps hj are given by (fT3j) . For j = 2 ^ gives 

<t>2,n(hl(u , ... , U k ), ■■■ , h m (u , ... , U k )) = h 2 (f n (u , ... , U k ),U , ... , U k -l) 

= no * h(ui,U2 ■■■ , u k ^ m+ \) 
= hi(u , ■ ■ ■ ,u k ). 

Therefore, <f>2 n projects into coordinate 1. Generally, for j > 2 we have 

<t>j,n(hi(u , . . . ,u k ), . . . ,h m (u , . . . ,U k )) = hj(f n (u , ■ ■ .,Uk),U , ■ ■ -,U k -l) 

= Uj-2 * h{Uj-\,Uj . . . ,Uj + k-m-\) 

= hj-x(u , ■ ■ ■ ,u k ). 
Therefore, for each n and for every (t\, . . . , t m ) G H(G k+1 ) we have 

$ n (tx, . . . , t m ) = [g n (ti, . . . , t m ),t\, . . . , t m -i] 
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i.e., $ n \H(G k+1 ) is of scalar type. Since by Part (b) H(G k+1 ) = G m it follows that is of scalar 
type. ■ 

We point out that the SC factorization in Theorem E^a) does not require the determination of 
<j)j^ n for j > 2. However, as seen in Parts (b) and (c) of the theorem the rest of the picture fits 
together properly. 

We discuss several applications of Theorem [6] in later sections below. The following example 
gives an application in finite settings. 

Example 7 Let G be any nontrivial abelian group (e.g., Z2 in Example viewed as an additive 
group) and consider the following difference equation of order two 

x n+ i = x n -i + a, a,x ,x-i € G, a / 0. (21) 

The unfolding of Eg. h21\) is the map F(u,v) = [v + a,u]. Note that the function H(u,v) = u + v 
is of type u + h(v) with h being the identity function on G (H is the same as Hq in Example 
With this H we have 

H(F(u, v)) = u + v + a = (f>(H(u, v)) 

where (ft : G — > G is defined as (ft(t) = t + a. Therefore, F > (ft. Using Theorem^ we obtain the SC 
factorization of Eq. \2l\i as 

tn+l = 4>(tn) =t n + a, t = H(x , X-l) = Xq + £_i 

3.3 Reduction types and chains 

Theorem [6] leads to a natural classification scheme for order reduction which we discuss in this 
section. We begin with a definition. 

Definition 8 The SC factorization of Eq.fcjfy into $16\) and Ji7| ) gives a type-(m, k + 1 — m) 
order reduction (or just type-(m,k + 1 — m) reduction) for ([!]). We also say that (OP is a 
type-(m, k + 1 — m) equation in this case. 

A second-order difference equation (k = 1) can have only the type-(l,l) order reduction into 
two first order equations. A third-order equation has two order reduction types, namely (1,2) and 
(2,1), a fourth order equation has three order reduction types (1,2), (2,2) and (2,1) and so on. Of 
the k possible order reduction types 

(l,fc), (2,fc-l),--. .(fc-1,2), (fc.l) 

for an equation of order k + 1 the two extreme ones, namely (1, k) and (k, 1) have the extra appeal 
of having an equation of order 1 as either a factor or a cofactor. 
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Eq.([T|) may admit repeated reductions of order through its factor equation, its cofactor equation 
or both as follows: 



Eq.© 



tn+l — Qnitni 



i tn- 



■m+1 



factor equation 

~ t n +l * h(x n , . . . , 2?n— Ai+m) 



factor — 
cofactor 

factor — 
cofactor 



cofactor equation 

In the above binary tree structure, we call each branch a reduction chain. If a reduction chain 
consists only of factor (or cofactor) equations then it is a factor (or cofactor) chain. In particular, 
we show later (Corollary 1201 below) that a linear difference equation of order k + 1 has a full cofactor 
chain leading to a system of k+1 linear first order equations, as discussed earlier in the introduction. 
The next example exhibits a full factor chain. 



Example 9 Consider the following reduction chain: 



X n +1 — %n ~\~ 



a(x n x n —i) 



t 



ati 



n+1 



Xn—1 Xn—2 



Sn+1 
t n +l 



as n 



(22) 



tn-l 

Xn+1 — tn+l ~i~ X n 

The 3rd order equation is reduced by substituting t n = x n — x n -\, then the 2nd order factor equation 
is reduced by the substitution s n = t n /t n -\ (see Definition [73l below and the comments that follow 
it). The factor chain in this example has length three as follows: 



Xn+1 — X n t 



a{x n x n — i) 



tn+l 



ati 



s n+l 



as r 



Xn—l X n —2 t n —i 

Since a first order difference equation does not have order-reducing form symmetries the SC 
factorization process of an equation of order k + 1 must stop in at most k steps. The result then 
is a system of k + 1 first order equations that mark the ends of factor/cofactor chains. In (|22[) the 
three first order equations that mark the ends of factor/cofactor chains can be arranged as follows: 



s n+l 

t n +i 

X n +1 



as n 

&Sntn 

&Sntn ~i~ X n . 



(23) 



Since each equation depends on the variables in the equations above it, (|23p is a "triangular 
system." In principle, a triangular system can be solved by solving the top-most equation and then 
using that solution to solve the next equation and so on. In the case of (|23p since the equations are 
linear (including the nonautonomous and nonhomogenous versions) an explicit formula for solutions 
of the system can be determined; in particular, the third order equation (|22p is integrable. Whether 
a complete SC factorization as a triangular system exists for each higher order difference equation of 
type (H|) is a difficult question that is equivalent to the apparently more basic problem of existence 
of semiconjugate relations for ([T|). 
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4 Reduction of order of difference equations 



In this section we use the methods of the previous section to obtain reductions in orders of various 
classes of difference equations. 

4.1 Equations with type- (A;, 1) reductions 

For type- (A;, 1) reductions of Eq.([T]) m = k. Therefore, the function h : G — » G in (| 13j) is of one 
variable and yields the form symmetry 

H(u , . . . , u k ) = [u * h(ui), u\ * h(u 2 ) u k _ x * h(u k )]. (24) 
Theorem [6] gives the SC factorization as the pair 

tn+l = 9n(tn,t n -l, . . . ,t n -k+l) (25) 
X n +1 = t n +l * [hiXn)]' 1 . (26) 

The functions g n : G k — > G are determined by the given functions /„ in as in the previ- 
ous sections. From the semiconjugate relation it follows that a type-(A, 1) reduction with a form 
symmetry of type (|24p exists if and only if there are functions g n such that 

g n (u * h(ui),m * h(u 2 ), ... , u k ~i * h(u k )) = f n (uo, u\,... ,u k ) * h(u ). (27) 

Remark 1. The order of the factor equation (|25p is one less than the order of (pQ). If a solution 
of (|25p is known then by the factorization theorem the corresponding solution x n of (pQ) is 
obtained by solving the cofactor equation (|26p . which has order one. For this reason we may refer 
to Eq. ([25j) as an order reduction for Eq.([T]). 

4.1.1 Invertibility criterion 

The next result from [19] gives a necessary and sufficient condition for the existence of <j) n when 
the function h above is invertible. We repeat the proof for the reader's convenience. 

Notation: the symbol /i -1 (-) denotes the inverse of h as a function and the symbol h(-) 
denotes group inversion. 

Theorem 10 (Invertibility criterion) Let h : G — > G be a bijection. For uq,vi, . . . ,v k € G let 
Co = uq and define 

Q = h- 1 (^\*v j ), j = l,...,k. (28) 

Then Eq.^j has the form symmetry \2J$ with SC factors (f> n satisfying Ji?7| ) if and only if the 
quantity 

fn{u ,(i,. ■ ■ ,(k) *h(u ) (29) 

is independent of uq for every n. 
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Proof. First assume that the quantity in f|29|) is independent of u$ for all v\,. . .,v k so that the 
function 

g n (y x , . . . ,v k ) = f n (u , Ci, ■ ■ ■ >Cfc) * h(u ) (30) 
is well defined. Next, if H is given by (|24|) then for all no, u\, . . . ,u k 

g n (H(uQ,ui,...,u k )) = g n (uQ*h(ui),ui * h(u 2 ), ■ ■ ■ , u k -i *h(u k )). 

Define 

Vj = Uj-i*h(uj), j = l,...,k. (31) 

Then by 

Ci = /i _1 (n " 1 * v\) = /i _1 (uq 1 * uq * h(ui)) = u\. 
In fact Q = Uj for every j, for if by way of induction Q = ui for 1 < I < j then 

Cj = h~ x {C^} x * vj) = h~ x {u'T^ 1 * Uj-x * h(uj)) = uj. 

Now by (pfflj) 



g n (H(u ,ui, . . .,u k )) = f n (u , ... ,u k )* h(u ) 
so if F n and $ n are the unfoldings of f n and g n respectively, then 

H(F n (u , . . .,u k )) = [f n (u , ... ,it fe ) * h(u ),u * h(ui), . . ■ ,u k _ 2 * h(ujfc-i)] 

= bn(#(«0, «1, • • • , «fc))> «0 * ■ • • , «fc-2 * ^(«Jfc-l)] 

= $ n (H(u ,...,u k )) 

Finally, to show that H is a semiconjugate link for Eq.([T]) we show that H is onto G k . Let 
(vi, . . . , v k ) E G k and let u k be any element of G, e.g., the identity. Define 

n i-i = u i * [^(""j)] -1 ) 3 = k, k - 1, . . . , 1 

and note that 

H(u , . . . , u k ) = [uq * /i(ni), mi * h(u 2 ) u k _i * h(u k )\ 

= [U1,U2, ■•• 

i.e., -ff is onto G fe as claimed. 

Conversely, if H as given by (|24p is a form symmetry then the semiconjugate relation implies 
that there are functions g n such that 

f n (uo, ■ ■ ■ , u k ) * h(uo) = g n (u * h(ui), u k -i * h(u k )). (32) 
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For every v \ , . . . , v k in G and with Q as defined in (|28|) , 

fn(u , Cl, • • • , Cfe) * UO = 9n{uo * /l(Cl), Cl * • • • , Cfc-1 * KCfc)) 

= g n {vi, ...,v k ) 

which is clearly independent of uq. ■ 

It is worth noting that (|28p is a backwards version of the cofactor equation (|26p that is obtained 
by solving it for x n instead of x n +i- To do this we required h to be invertible. Of course, in ()28[) it 
is necessary to iterate only k times. 

4.1.2 Identity and inversion form symmetries 

We now examine two of the simplest possible form symmetries within the context of the previous 
section that are based on the scalar maps 

h\(u) = u and = u~ . 

Each of h±\ is invertible and equals itself (self-inverse maps). 

Definition 11 The form symmetry of type \2J$ that is generated by hi is the identity form 
symmetry and that generated by is the inversion form symmetry. 

The next result is an immediate consequence of Theorem [TU1 

Corollary 12 (a) For every uq, v±, . . . , v k G G let Qq = uq and define Cj = Cj\ * v j for j = 1, . . . , k. 
Then Eq. FT]) has the identity form symmetry if and only if the quantity 

fn(uo,(i,---,(k)*u (33) 

is independent of uo for every n. 

(b) For every uo,v\, . . . ,Vk & G let £o = an d define (j = vj 1 * (j-i for j = 1, . . . , k. Then 
Eq. {TJ) has the inversion form symmetry if and only if the quantity 

fn(u ,Ci,---,Ck)*u 1 (34) 

is independent of uq for every n. 

For the inversion form symmetry we can improve on Corollary|12lfb) considerably. In fact, there 
is a simple characterization of functions f n which satisfy (|34j) . 

Definition 13 Equation |7]) is homogeneous of degree one (HD1 ) if for every n = 0, 1, 2, . . . 

the functions f n are all homogeneous of degree one relative to the group G, i.e., 

f n (u * t, . . . , U k * t) = f n (u , ... ,U k )*t 

for all t,U{ £ G, i = 0, . . . , k, and all n > 0. 
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Note that the third order equation in Example [9] is HD1 with respect to the additive group 
of real numbers while its second order factor is HD1 with respect to the multiplicative group of 
nonzero real numbers. 

There is an abundance of HD1 functions on groups. For instance, if G is a nontrivial group, k 
a positive integer and g : G k — > G is any given function, then it is easy to verify that the mappings 
gj : G k+1 -» G defined by 

gj(u , ...,u k )= g{u * u^ x ,u\ * u^ 1 , . . . ,u k -i * u^ 1 ) * uj 

are HD1 functions for every j = 0, 1, . . . , k. Further, if gi,g 2 ■ G k+1 — > G and / : G 2 — > G are HD1 
then so is the composition 

f(gi(u , ■ ■ .,u k ),g 2 (uo, • • -,u k )). 

The next result was proved by direct arguments in |18j : also see [22]. Here we give a different 
proof using Corollary [T2Tb) and Theorem [6J 

Corollary 14 Eq. (Tip has the inversion form symmetry if and only if f n is HD1 relative to G for 
all n. In this case, |7]j has a type-(k,l) order-reduction with the SC factorization 

tn+l = /n(l,*n > (*n * *n-l)~\ • • • , (*n * *n-l * ■ • • * in-fe+l) _1 ) (35) 
£n+l = t n +l * (36) 

Proof. First if Q = vj 1 * Cj-i as i n Corollary I12lfb) then by straightforward iteration 

Cj = vj 1 *■■■ * v^ 1 *u , j = 1, . . . , m. (37) 

Now if f n is HD1 for every n then 

fn(uo,(u ■ ■ ■ Xk) * Uq 1 = / n (l,Ci * Uq 1 , . . . ,(k * Uq 1 ) 

which by (|37p is independent of uq. Thus by Corollary|12lfb) Eq.([I]) has the inversion form symmetry. 

Conversely, assume that ([1]) has the inversion form symmetry. Then by Corollary [T2T b) for 
every Uq,v±, . . . ,v k G G the quantity in ([M]) is independent of no- Thus there are functions <p n 
where 

fn(u ,(i, ■■■Xk)*u 1 = <j) n (vi,.. .,v k ) (38) 

= </>„( Co * Cf" 1 , Cfc-i ^aT 1 )- 

Note that (j38j) holds for arbitrary values of uq, Cl> • • • > Cfc since «i, . . . , are arbitrary. Thus for 
all t, so, • • • , s k S G and all n, 

/ n (s * t, . . . , S k * t) = 4> n {(so *t)* (Si * ty 1 , (s k -l *t)* (s k * ty 1 ) * (s * t) 

= [(f> n (s * sj 1 , S k -l * s^) * s ] * t 

= fn(so, ■ ■ ■ ,S k ) *t. 
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It follows that f n is HD1 relative to G for all n. 
Finally, the inversion form symmetry 

H(uq, ...,u k ) = [u * u^ 1 , . . . ,u fc _i * u^ 1 ] 

yields a semiconjugate relation that changes variables to t n = x n /x n _\ and the SC factorization 
([35]) and ([36]) is obtained by Theorem [6] ■ 



Example 15 Consider the 3rd order difference equation 

x n +i = —t > a n ,b n > 0. (39) 

Here for all n, f n (uo,ui,U2) = uqU\/ (a n UQ + b n U2) is HD1 on (0, oo) under ordinary multiplication. 
By Corollary [77] Eq. \39\) has the following SC factorization 

""nJ"nJ"n—X ~+~ "n 
%n+l — t n +\X n . 

Next, consider the factor equation fl^0| ) again on (0, oo) under ordinary multiplication. This equation 
has order two and represents the order reduced form of Eq. t fgpj) in the sense of Remark 1. Note that 
(70[ j is not HD1; however, substituting £i = vi/uq in the functions <f> n (uo,ui) = u\/(a n UQUi + b n ) 
gives 

( V\\ vx/uq Vi 

<Pn [UO, ) U = ; ■ — 7—Uq - 



uq) a n u vi/u + b n a n vi+b n 

which is independent ofuQ. Thus Corollary [T^ implies that Eq. fl^0|) has the identity form symmetry. 
Its SC factorization is easily found using Theorem^ upon substituting r n = i n i n -i to get 

r n+ i = \ h (41) 
. _ r n+l 

tn+l — — • 

The further substitution s n = l/r n in (JW produces a linear equation s n+ i = a n + b n s n and yields 
a full factorization of Eq. ( f35j) into the following triangular system of first order equations 

Sn+l = On "I" b n S n 
1 



t n +l ~ 
%n+l — t n -^\X n . 
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The factorizations in the preceding example were used in analyzing the global behavior of 
Eq. ([39|) in [20]. SC factorization of HD1 equations has also been used in [5] and [12] to study the 
global behavior of difference equations of order two. 

Corollary 16 Let T be a nontrivial algebraic field and a\, . . . ,a k € T with a k ^ 0. Each linear 
homogeneous difference equation 

%n-\-l — &l%n + ' ' ' + dk^n—k 

of order k + 1 is HD1 on the mulitplicative group .F\{0}. Thus it admits a type-(k, 1) reduction 
with factor equation 

tn+l =0,1 + - 1 h — " • 

f-n ^"nXn—1 ' ' ' ^n— fc+1 

The factor equation in Corollary [16] is known as the Riccati difference equation of order k; see 

a- 

4.1.3 The linear form symmetry 

In practice, the one-variable function h is often defined using structures that are more complex 
than a group. In particular, on a field T the function 

h(u) = —au 

where a is a fixed nonzero element of the field defines a form symmetry 

H(u ,ui, . . . ,u k ) = [uq - au\,u\ - au 2 , . . . ,u fc _i - au k ]. (42) 

For convenience, we represent the field operations here by the ordinary addition and multiplica- 
tion symbols. We call H defined in (j32"|) the linear form symmetry. Relative to the additive group 
of a field, the linear form symmetry generalizes both the identity form symmetry (a = —1) and the 
inversion form symmetry (a = 1). 

The linear form symmetry is characterized by a change of variables to 

The factor and cofactor equations for the SC factorization of ([T]) with the linear form symmetry 
are given by (|25|) and (f26l) as 

t n +l = ^n(*n> • • • j*n-fc+l)> (43) 
Xn+l = t n +l + OiX n (44) 
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For a given sequence {t n } in T the general solution of the linear cofactor equation (|44p is easily 
found to be 

n 

x n = a n x + ^2 ® n ~ j tj, n > 1. (45) 

3=1 

This equation gives a solution of Eq.([T]) if {t n } is a solution of (|43p with the form symmetry 
(|42p . Difference equations possessing the linear form symmetry have been studied in [6] and [21] 
where the SC factorization above has been used to reduce equations of order 2 to pairs of equations 
of order 1. 

The following corollary of Theorem [10] gives a condition for verifying whether Eq.([T]) has the lin- 
ear form symmetry. For easier reading we denote the multiplicative field inversion by the reciprocal 
notation. 

Corollary 17 For arbitrary uq, v%, . . . , u& in a field T define Co = uq and 

i=l 

Then Eq. |ip /ias i/ie linear form symmetry (f^SD i/ and only i/ i/ie quantity 

fn(uo,(i, - ■ ■ ,(k) -au 

is independent of uq. 

Note that (|46p defines in Corollary 1171 explicitly rather than recursively. It is obtained from 
the recursive definition (|28p by a simple calculation; since 



equality ()46p can be established by direct iteration. The next example illustrates Corollary 1171 

Example 18 Let J 7 be a field and g n : T — > T be a given sequence of functions. Consider the 
difference equation 

x n+ i = ax n -j + g n (x n + cx n _ k ) (47) 
where c^0,0<j<k and we define f n (uo, u±, . . . , u^) = auj + g n (uQ + cuk)- Using ffljfy we obtain 
/ k \ 



i ( 
a \ v-^ avi I r c i ^— «, 



i=l V i=l 



a fc-i+l 
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If there is a £ J 7 such that 

-a + -^r = or a = a j+1 and c = -a k (48) 
then by Corollary [7?] Eq. |^7| ) has the linear form symmetry with a SC factorization 

tn+l = ~ J] ait n-i+l +9n\Y^ O^n-f+l ) , ( 4 9) 
i=l \i=l / 

n 

x n+ i = t n+ i + ax n = a n x + ^ a n ~ j tj. 

3=1 

Remark. The above example in particular shows that the nature of the field T is important in 
obtaining SC factorizations. For example, the difference equation 

satisfies (1351) on Z2 but not on C. 



4.2 Equations with type-(l, k) reductions 

In this section we obtain type-(l, k) reductions for a class of equations on the field C of complex 
numbers that includes all linear nonhomogeneous difference equations with constant coefficients as 
well as some interesting nonlinear equations. 

For a type-(l, k) reduction m = 1. Therefore, h : G k — > G and the form symmetry has the scalar 
form 

H(u , . . . , u k ) = n * h{u\, u k ). (50) 
This form symmetry gives the SC factorization 

tn+l = 9n(tn) (51) 

x n+ i = t n+ i * [h(x n , . . . , x n _ fc+ i)] _1 . (52) 



Remark 2. The order of the cofactor equation (|52p is one less than the order of {[TJ) . We can obtain 
the solution £ n+ i of the factor equation (|25() . which has order one, either explicitly or determine its 
asymptotic properties. For this reason we may refer to Eq. ([52|) as an order reduction for Eq.([T]). 

4.2.1 Form symmetries and SC factorizations of separable equations 

Consider the additively separable difference equation 

Xn+l = a n + 4>s{Xn) + 4>l(Xn-l) H 4>k{x n -k) (53) 
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where k > 1 is a fixed integer and (|53p is defined over a nontrivial sub field T of C, the set of 
complex numbers. Specifically, we assume that 

i-j,a n ef, 0j : T -> J 7 , j = 0, 1, ... n = 0,1,2, .. . (54) 

We use subfields of C rather than its subgroups because the results below require closure under 
multiplication. The subfield that is most often encountered in this context is the field of all real 
numbers R. The unfolding of Eq. (|53p is defined as 

F n (z ,z 1 ,...,z k ) = [a n + (f)o(zo) H \-(j) k (zk),za,...,z k „x]. 

Based on the additive nature of Eq. (|53j) we look for a form symmetry of type 

H(z Q , zi, . . . , z k ) = Zq + /ii(zi) H h h k (z k ) (55) 

where hj : T — > C for all j; we do not require that T be invariant under the maps hj. 
With these assumptions, equality © takes the form 

a n + fa(z ) H 1- (Pk(z k ) + ^i(^o) H 1- h k (z k -i) = g n (zo + /ii(^i) H h h k (z k )) (56) 

where : C — > C. Note that if (|56p holds and the subfield T is invariant under hj for all j then T 
is invariant under g n also; i.e., 5 f n (-^ r ) C F for all n. 

Our aim is to determine the functions hj and g n that satisfy the functional equation (I56j) . To 
simplify calculations first we assume that the functions fa and hj are all differentiable on T. The 
differentiability assumption can be dropped once the form symmetry is calculated. 

Take partial derivatives of both sides of (|56p for each of the k + 1 variables Zq, zi, . . . , z k to 
obtain the following system of k + 1 partial differential equations: 

</>d(zo) + h[( z o) = g'nizo + hi(zi) -\ h /ifc(zfc)) (57) 

<t>'i{zx) + h' 2 (zi) = g' n {zo + hi(zi) + ■■■ + h k (z k ))h[(zi) (58) 

0fc-i(^fc-i) + ^fc(^fc-i) = 9n( z o + /ii(^i) H h M^jO^fe-iC^-i) 

^jii(^) = 0n(*O + ^lOl) H 1- h k (z k ))h' k (z k ). (59) 

Taking equations ([57]) and ([5"8]) together, we may eliminate 5^ to obtain 

4Wo) + fti(«o) = 777 \ • (60) 

The variables zo and z\ are separated on different sides of PDE (I60p so each side must be equal 
to a constant c G C. The left hand side of ([60]) gives 

fa (z) + = c h\{z) = cz - 4>a{z). (61) 
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Similarly, from the right hand side of (|60p we get 

(friiz) + h' 2 (z) 

' = c h 2 (z) = chi(z) - fa(z). 

Note that from ([57]) and (JHU) it follows that 

g' n {z + h 1 (z l ) + --- + h k {z k ))=c. (62) 

It is evident that this calculation can be repeated for equations 3 to k — 1 in the above system 
of PDE's and the following functional recursion is established by induction 

hj = chj-i- 4>j-i, j = 1, . . . ,k, h (z) = z. (63) 

Finally, the last PDE ([59]) gives 

ch k = <]) k . (64) 

On the other hand, setting j = k in ([63]) gives h k = ch k _i — 4> k ~i so equality ([Ml) implies that 
there must be a restriction on the functions 4>j. To determine this restriction precisely, first notice 
that by ([MD and ([63]) 

c 2 h k ^i(z) - c0 fc _i(z) - 4> k (z) = for all z. 

Now, applying (|63|) repeatedly k — 1 more times removes the functions hj to give the following 
identity 

c k+1 z-c k M z ) -c k ~ l (j)i{z) c0 fc _i(z) -<t> k {z) =0 for all z. (65) 

This equality shows that the existence of form symmetries of type (|55p requires that the given 
maps (f>o, 0i, . . . , (f>h plus the identity function form a linearly dependent set. Equivalently, ([65]) 
determines any one of the functions 4>q , </>i , . . . , 4> k as a linear combination of the identity function 
and the other k functions; e.g., 

<j) k {z) = c k+1 z - c k Mz) - c^M*) c0 fc _ a (*). (66) 

Thus if there is c S C such that (|65p or (|66p hold then Eq. (|53p has a form symmetry of type 
(|55p . This form symmetry H can now be calculated as follows: Using the recursion (|63p repeatedly 
gives 

hj(z) = c*z - c*~ x 4)q(z) (pj-i(z), j = l,...k. (67) 

Therefore, 

H(z ,zi, ...,z k ) = z + h^zt) H h h k (z k ) 

k 

= Zo + ^2,[c j Zj - C j ' l (t)Q{Zj) <£,-_i(Zj)]. 
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With the form symmetry calculated, we now proceed with the SC factorization of Eq. (|53p . To 
determine the first order factor equation, from (|62p we gather that g n (z) is a linear non-homogeneous 
function of z. The precise formula for g n can be determined from (|56|) . (|63p and (|64|) as follows: 
First, note that if the constant c is in the subfield J- then by (|67p T is invariant under all hj and 
thus also under all g n . In this case, the function H(zq, z±, . . . ,-Zfc) is onto JF, i.e., i?(.F fc+1 ) = JF. 
This is easy to see, since for each z G T and arbitrary zi,...,z^ G J- we may choose zo = 
2 — [h\{zi) + • • • + hk(zk)] £ J 7 so that z = H{zq, Z\ y . . . , z k ). On the other hand, if c £ T then a 
similar argument proves that H is onto C. Now, we calculate g n on T or on all of C as 

9n(z) = 9n(z + /ii(^i) H h h k (z k )) 

k 

= a n + ^[^'-i(^'-i) + hj(zj-.!)] + <f>k(zk) 
0=1 

= a n + c[z + /ii(zi) H h /ifc(zfc)] 

= a n + cz. 

This determines the factor equation (]5ip for Eq. (l53p : the cofactor equation is given by (I52p . 
Therefore, the SC factorization is 

z n+1 = a n + cz n , z = xq + h\{x-i) H h h k (x_ k ) (68) 

a^n+i = ^n+i - hi(x n ) h k {x n _ k+ i). (69) 

Note that the functions hj do not involve the derivatives of an d we showed that the form 
symmetry H is a semiconjugate link without using derivatives. Therefore, the SC factorization is 
valid even if <pi are not differentiate. The following result summarizes the discussion above. 

Theorem 19 Eg. t53\) has a separable form symmetry 155\) if there is c G C such that \65\) holds. 
In this case, 153\) reduces to the system of equations 168\) and 169\) of lower orders with functions 
hj given by (67| ). 

4.2.2 Complete factorization of linear equations 

A significant feature of Eq. ([69]) , which by Remark 2 represents a reduction of order for Eq. (|53p , is 
that ([69]) is of the same separable type as ([53]) . Thus if the functions hi, . . . , h k satisfy the analog 
of (I65p for some constant d G C then Theorem 1191 can be applied to (I69p . In the important case 
of linear difference equations this process can be continued until we are left with a system of first 
order linear equations, as seen in the next result (for its proof and related comments see |22j) 

Corollary 20 The linear non-homogeneous difference equation of order k + 1 with constant coef- 
ficients 

x n+1 + b x n + &ix n _i H h b k x n -i = a n (70) 
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where &q, o n 6 C is equivalent to the following system of k + 1 /irsi order linear non- 

homogeneous equations 



Z0,n+1 — OL n + C()Zo,ni 
%l,n+l = Zo,n+l + c l z l,n 



Zk,n+1 — %k— l,n+l "I - Ck%k,n 




The preceding result shows that the classical reduction of order technique through operator 
factorization is also a consequence of semiconjugacy. Indeed, even the eigenvalues can be explained 
via semiconjugacy. Corollary [20] also shows that the complete factorization of the linear equation 
(|70p results in a cofactor chain of length k + 1. 

4.2.3 Multiplicative forms 

As another application of Theorem [19] we consider the following difference equation on the positive 
real line 



Pn,y-j G (0,oo), ipj : (0, oo) -> (0,oo), j = 0, . . . fc. 

This difference equation is separable in a multiplicative sense, i.e., with (0, oo) as a multiplicative 
group. Taking the logarthim of Eq.(|72|) changes it to an additive equation. Specifically, by defining 

x n = lny n , y n = e Xn , ln/3 n = a n , (pj(r) = lmpj(e r ), j = 0, . . . k, r£l 

we can transform (|72p into (|53[) . Then Theorem 1191 implies the following: 

Corollary 21 Eq.^7^) has a form symmetry 




Vn+l = Mo(y n )lpl{yn-l) ■ ■ ■ IpkiVn-k), 



(72) 




H(to,h 



. . . ,tk) = tohi(t\) ■ ■ ■ hk{tk) 



(73) 



if there is c £ C such that the following is true for all t > 



MtfMt) 



■■■Mt) = t' 



c k+i 



(74) 
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The functions hj in 



7 Hfy are given as 



hj(t) = i c > (i)- cJ_i • • • ^i-i(t), 3 = 1, • • • k (75) 



and the form symmetry in | 75) and |75[) yields the type-(l, k) order reduction 



r n+ i=(3 n r c n , r = y Q hi(y-!) ■ ■ ■ hk{y-k) (76) 

2/n+l = T— ( \ 7—7 s ■ (77) 

hi{y n ) ■ ■ ■ hk{yn-k+i) 

For a different and direct proof of the above corollary when k = 1 see |23| . The next example 
applies Corollary [21] to an interesting difference equation. In particular, it shows that the global 
behavior of a population model discussed in [8] is much more varied than previously thought. 

Example 22 (A simple equation with complicated multistable solutions). Equations of 
type 153\) or are capable of exhibiting complex behavior, e.g., coexisting stable solutions of many 
different types that range from periodic to chaotic. Consider the following second-order equation 

x n+ i = x n - X e a - Xn - Xn -\ x-uxq > 0. (78) 

Note that Eq. fl7<§l ) has up to two isolated fixed points. One is the origin which is repelling if a > 
(eigenvalues of linearization are ±e a / 2 ) and the other fixed point is x = a/2. If a > 4 then x is 
unstable and non-hyperbolic because the eigenvalues of the linearization of (f7ffp are —1 and 1 — a/2. 
Figure 1 shows the bifurcation of numerically generated solutions of this equation with a = 4.6. 
One initial value x-i = 2.3 is fixed and the other initial value xq ranges from 2.3 to 4-8; i.e., 
approaching (or moving away from) the fixed point x = 2.3 on a straight line segment in the plane. 

Stable solutions with periods 2, 4, 8, 12 and 16 can be easily identified in Figure 1. All of 
the solutions that appear in this figure represent coexisting stable orbits of Eq. (75|). There are also 
periodic and non-periodic solutions which do not appear in Figure 1 because they are unstable (e.g., 
the fixed point x = 2.3J. Additional bifurcations of both periodic and non-periodic solutions occur 
outside the range 2.3-4-8 which are not shown. 

Understanding the global behaviors of solutions of Eq. {78fy is relatively easy if we examine its 
SC factorization given by (761) and (770. Here k = 1 and 



ipo(t) = e *, ipi(t) = te *, p n = e a for all n. 

Thus ( f7^[ ) takes the form 

ip (t) c i/ji(t) = t° 2 for allt>0 
e^He' 1 = t c2 for all t > 
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Figure 1: Bifurcations of solutions of Eq. ([78|) with a changing initial value; a = 4.6 is fixed. 
The last equality is true if c = — 1, which leads to the form symmetry 

-1 , u\-l 1 . rr, .. \ U 



hi{t) = t~Vo(i) = —— => H(uq,ui) 



te t u±e Ul 

and SC factorization 



r n +i = —, r = x hi(x-i) = X ° (79) 

r n x-ie x - 1 

= TT^ = r n +ix n e- Xn . (80) 



All positive solutions of (791) with ro ^ e a l 2 are periodic with period 2: 



e" I I xq X-\e a x ~ 1 



r'o J I x^ie x ~ 1 xq 



Hence the orbit of each nontrivial solution {x n } of |7ff| ) in the plane is restricted to the pair of 
curves 

e a 

= — te"* and £ 2 (t) = rote' 1 . (81) 
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Figure 2: Two of the orbits in Figure 1 shown here on their loci of two curves £1, £2 in the state space. 

Now, ifx-i is fixed and xq changes, then tq changes proportionately to xq. These changes in initial 
values are reflected as changes in parameters in [80]). The orbits of the one dimensional map bte~ l 
where b = r$ or e a /r$ exhibit a variety of behaviors as the parameter b changes according to well- 
known rules such as the fundamental ordering of cycles and the occurrence of chaotic behavior 
with the appearance of period- 3 orbits when b is large enough. Eg. h80\) splits the orbits evenly over 
the pair of curves \81\) as the initial value xq changes and coexisting, qualitatively different stable 
solutions appear. 

Figure 2 shows the orbits of {78fy for two different initial values xq with a = 4.6; the first 100 
points of each orbit are discarded in these images so as to highlight the asymptotic behavior of each 
orbit. 

5 Concluding remarks 

In this article we covered the basic theory of reduction of order in higher order difference equations 
by semiconjugate factorization. Needless to say, a great deal more remains to be done in this area 
before a sufficient level of maturity is reached. It is necessary to generalize and extend various ideas 
discussed above. In addition, a more detailed theory is necessary for a full discussion of type-(m, k) 
reductions for all m. 

The list of things needing proper attention is clearly too long for a research article. An upcoming 
book [17J is devoted to the topic of semiconjugate factorizations and reduction of order. This book 
consolidates the existing literature on the subject and offers in-depth discussions, more examples, 
additional theory and lists of problems for further study. 
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